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Abstract 

We show that a finite Josephson Junction (JJ) chain, ending with two bulk super- 
conductors, and with a weak link at its center, may be regarded as a condensed 
matter realization of a two-boundary Sine-Gordon model. Computing the partition 
function yields a remarkable analytic expression for the DC Josephson current as a 
function of the phase difference across the chain. We show that, in a suitable range 
of the chain parameters, there is a crossover of the DC Josephson current from a 
sinusoidal to a sawtooth behavior, which signals a transition from a regime where 
the boundary term is an irrelevant operator to a regime where it becomes relevant. 
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1 Introduction 



In this paper, we analyze a superconducting 1+1-dimensional system, defined 
on a finite interval of length L. If the bulk is described by a massless the- 
ory, and conformal boundary conditions are chosen, one could understand the 
properties of the model, using the formalism of boundary conformal field the- 
ories [1]. If one deviates from this situation by either adding an interaction 
in the bulk, or at the boundary, or both, the behavior of the system becomes 
much more interesting, since it involves crossovers depending on the bulk and 
boundary energy scales, as well as on the finite size L. 
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In the sequel, we shall show that a JJ-chain with a weak link at its center 
and ending with two bulk superconductors at fixed phase difference (/?, is the 
prototype of a condensed matter realization of a two-boundary Sine-Gordon 
model [2], whose Hamiltonian is given by 
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Boundary field theories appear to be relevant in several different contexts. In 
condensed matter physics, they are mostly generalizations of quantum impu- 
rity models, which may be described by using the Tomonaga-Luttinger Liquid 
(TLL)-paradigm [3]; for instance, boundary interactions appear in the anal- 
ysis of the Kondo problem [4], in the study of a one-dimensional conductor 
in presence of an impurity [5] , and in the derivation of the tunneling between 
edge states of a Hall bar [6]. The TLL paradigm shows that many interactions 
are simply diagonalizable in the basis of appropriate collective bosonic modes, 
and that non diagonalizable interaction usually correspond to exactly solvable 
Hamiltonians, such as Sine-Gordon models [7,8]. Recently, boundary field the- 
ories have been investigated in the context of string theories. For instance, in 
studying tachyon instabilities [9], one is faced with the fact that the space 
of interacting string theory [10] may be mapped onto the space of boundary 
perturbations of conformal theories [11], and that the renormalization group 
fiow determined by boundary perturbations may be identified with tachyon 
condensation [12]. Affleck and Ludwig [13] showed that the boundary entropy 
g is decreasing along the renormalization group trajectories, triggered by the 
boundary interaction. 

In a inspiring paper [14], Glazman and Larkin analyzed the quantum phase 
diagram of a JJ-chain in the Vg — J-pIane, where Vg is an external gate voltage 
applied to each junction, while J is the Josephson coupling between neighbor- 
ing grains. They found evidence that this system undergoes a phase transition 
between an attractive TLL phase, with g < 1, and a repulsive TLL phase, 
with g > 1. While the former phase is the one-dimensional analog of the 
superconducting phase [15], the repulsive Tomonaga-Luttinger phase is pe- 
culiar of a one-dimensional system [3]. To be self-contained, here we shall 
provide a detailed field-theoretical description of the one-dimensional infinite 
chain analyzed in Ref.[14]: our rather pedagogical derivation evidences how 
the one-dimensional JJ-chain may be described in terms of interacting 1-1-1- 
dimensional chiral fermions and how, using the TLL paradigm, the interaction 
is exactly diagonalized in a pertinent basis. The TLL-(7 parameter [16] is cru- 
cial for the analysis of the phase diagram. Indeed, while for < 1 the system 
supports an attractive TLL phase (superconducting), ioi g > 1 the JJ-chain 
is described by a repulsive TLL phase [14]. The g — 1-line corresponds to 
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a noninteracting TLL model. All this features may be quantitavely derived 
within the framework of the bosonized 1+1-dimensional TLL-modeP . 

Using the TLL paradigm, we show that a finite JJ-chain with a weak link at its 
center is mapped onto a two-boundary Sine-Gordon model, with fixed Dirichlet 
boundary conditions at the outer boundary, and with dynamical boundary 
conditions at the inner boundary. To study the effects of the interaction at 
the inner boundary, we perform a Renormalization Group (RG) analysis, to 
derive how the effective parameters of the system scale as a function of the 
size of the chain L. We find that in the repulsive TLL- phase {g > 1) the 
boundary term is perturbative for any size L. At variance, in the attractive 
TLL-phase {g < 1), wc find that there is an RG-invariant length, L*, such 
that, for L < L* the boundary term is perturbative, while for L > L* it 
becomes nonperturbative. As for the models analyzed in [2], the crossover 
from the perturbative to the nonperturbative regime is evidenced by a change 
of the DC Josephson current (as a function of the phase difference at the bulk 
superconductors </?) from a sinusoidal to a sawtooth behavior. 

The paper is organized as follows: 

• In Section II we analyze the infinite one-dimensional JJ-array described in 
Ref . [14] and provide a detailed derivation of the mapping of this chain onto 
the anisotropic (XXZ) spin 1/2 model; 

• In Section III we construct the effective field theory for the equivalent XXZ 
chain. We bosonize the theory and identify the various parameters of the 
continuum model in terms of the microscopic parameters of the lattice 
Hamiltonian; 

• In Section IV, using the TLL-paradigm, we derive the phase diagram of the 
JJ-chain; 

• In Section V we show that the effective field theory for the JJ-chain with 
a weak link and ending with two bulk superconductors is indeed the two- 
boundary Sine-Gordon model; 

• In Section VI, using the Goulomb Gas Renormalization Group approach, 
we provide a careful estimate of the partition function of the two-boundary 
Sine-Gordon model for any value of g. We then derive the DC Josephson 
current across the chain, as a function of (p, at both fixed points and ex- 
plicitly show the existence of a crossover from a sinusoidal to a sawtooth 
behavior; 

• Section VII is devoted to a discussion of our results. 



Notice that here the TLL-g parameter is the inverse of the parameter used in 
Ref. [14]. 
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2 Mapping of the one- dimensional JJ-chain onto the XX Z spin-1/2 
model 



The simplest model Hamiltonian describing a one-dimensional JJ-chain is 
given by: 

E ^^'^ I 3 A/"\ ^ ^^'^ 
H^Hc + Hj^^Y.[-^9^^-Y) -•^Ecos(0,-0,+i) . (2) 



In Eq.(2) is the operator representing the number of Cooper pairs at 

site j in the phase representation and, thus, it takes only integer eigenvalues, 
rij, Ec is the charging energy of a grain, J is the Josephson couphng energy 
and N accounts for the influence of a gate voltage, since eN oc Vg. The sum 
over j ranges over the (L/a) sites, with L being the length of the chain, and a 
being the intergrain distance; imposing periodic boundary conditions amounts 
to fix (pL/a+j = 4>j- For J/Ec 0, the chain is an insulator for almost any 
value of N, since it costs an energy ~ Ec, to change the number of pairs at 
any grain: Ec measures, then, the insulating gap. 

When N = n + 1/2 + h, with integer n and /i ^ 1, the two states rij = n 
and Uj = n+1 become almost degenerate in energy, even for large Ec] in this 
hmit one may restrict the set of physical states to the Fock space F spanned 
by the 2^/« states 

L/a 

|{^}) = n l^j) ; rij ^ n,n + 1 . 

3=1 



The Josephson coupling hfts the degeneracy between n and n + 1, since Hj 
may be represented as 

Hj = --[e^'^i+ie"*'^^ + e''^^e-''^^+'] , (3) 
2 



with the operator e**^^' (e "^^) raising (lowering) the charge rij by 1. 

Resorting to the a well known procedure [17], one may easily construct the 
effective Hamiltonian Hes, describing the JJ-chain on the reduced space F [14]. 
Let P be the projector onto F and P± be the projector onto the subspace F±, 
to 0{J^/Ec), H^s takes the form: 



H,s^P{Hj + Hc)P + P 



HjP±Hj 

-JqEc 



p 



(4) 
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When restricted to F, the operators e^'^'^^ and —i^ niay be represented with 
the spin- 1/2 operators Sf and S^, as 

Pe^^^^P^Sf ;pLi^-n-i)p = 5| . (5) 



^From Eq.(5), one immediately sees that a charge n (n + l)-state corresponds 
to a spin-(-l/2) (1/2)- eigenstate of S^, and that, to 0{J'^/Ec), H^^s is given 

by 



To account for the contributions coming from intergrain capacitances, it is 
sufficient to retain only the nearest-neighbor terms [14] , since next-to-nearest 
neighbor hopping terms would give rise to irrelevant operators, and to add to 
the Hamiltonian (6) the term Ez Ejii S-S^^^ {Ez > 0) [14]. Thus, the system 
is usefully described [14] by 

He,-~ US^Si^i + - E ^- + A E S^S]^, , (7) 

^ j=l j=l j=l 



with H = Ech and A ^ Ez - ^i^. 

Eq.(7) is the Hamiltonian for a spin- 1/2 X X Z-chain in an external magnetic 
field H. The anisotropy parameter A may take positive, as well as negative 
values, depending on the constructive parameters of the JJ-chain. As eluci- 
dated in the following Sections, the sign of A is crucial for the emergence of a 
Repulsive Tomonaga-Luttinger (RTL) phase in a JJ-chain. 



3 Continuum field theory and bosonization of the XXZ-chain 



Using the standard bosonization technique [16], we derive in this Section the 
effective low-energy long-wavelength field theory associated to the Hamilto- 
nian (7). For this purpose, one starts to write the spin operators Sj in terms 
of Jordan- Wigncr (JW) spinless lattice fermions aj [18], obeying standard 
anticommutation relations: 

{aj,al}^5jk ■ (8) 
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The JW transformation amounts to define: 

i-i I i-i 

= a] exp[i7r ^ ajai] ; Sj = a]aj - - ; = exp[-i7r ^ alai]aj , (9) 
1=1 ^ i=\ 



which, in turn, imphes: 

[5';,5^] = z5,,e'^''^5'; . (10) 

^From Eqs.(9), the fermionic efi^ective Hamiltonian may be written as: 

LI a 

J L/a L/a L/a 

j=i j=i j=i 



The hopping term Hk is readily diagonahzed by resorting to the Fourier com- 
ponents of Qj, Gk, 



1 



a 



•J 



Y.^ke'"^"^ , (12) 



^JL/a k 



leading to: 

Hk = 'Ye{k)a\.ak ; (e(A;) = —Jcos{ka) — H) (13) 



If I if I < J, the Fermi surface is disconnected and consists of two isolated 
points at i/cF, with kp — arccos(if/ J). Keeping only the excitations about 
the Fermi points with momenta k such that |A; ± kp] < A, one obtains: 

|fe-feF|<A |fe+feF|<A |p|<A 

- Jsin(aA;F) ^ sin(pa)a|^(p)aR(p) , (14) 

|p|<A 



with: 

ttLip) = ttp+kp ; aR{p) = ttp-kp (bl < A) ■ 
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For A <^ kp, one may define the continuum chiral fields iPl/r{x) as 

-^^e''^^^^Mxj) + e-"'^"'Mxj) , (15) 



with Xj — ja; one gets, then 

Hk^ J sin(/cFa) ^ pWL{p)aL{p) - a^(p)aR(p)] 

bl<A 

= -^v/jd. [Ai^)"^ - ^k-)^] (16) 



where the Fermi velocity is given by vp = 2TrJ sm{akp). 

Eq.(16) is, of course, the effective low-energy theory of the hopping Hamilto- 
nian Hk', the cutoff A will be specified later. 

The dynamics of the fermionic fields ipL and ipu in the Heisenberg represen- 
tation, is described by 

M^, t) = + VFt) = ^ E , (17) 



and the equal time anticommutation relations are given by 

{Mp):i^l{p)} = ^,r>' ; {Mp)M)} = ; {Mp),i^lip')} = o -(is) 



Since ip^ip) with p > creates positive energy left-handed states, while ip^ip) 
creates positive energy right-handed states if p < 0, the "Fermi sea" fermionic 
ground state is defined as 

|FS) = ll[^l{p)ijU-p)m (MpM = MpM = 0) . (19) 

p<0 



Thus, by choosing A = l/(4a), one gets: 

{FS\2na[e-''''^'^^ijl{x,)M^j) + e''"'-^^ipU^j)M^jWS) = (20) 
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and: 

{¥^\a[^l^i{x,)M^,) + ^i(a;,)V^«(a;,)]|FS) = \ ; (21) 

thus, Sj is normal ordered respect to |FS), i.e., 

SI = 27ra[: ipl{xj)ipL{xj) : + : ipR{^j)ipR{xj) :] + 

2na[: ^l;l(xj)M^j) ■ e-'''-^^+ : : e'''-^^] (22) 

where :: denotes normal ordering. 

Using fermionic coordinates, one should now evaluate the Ising-Neel interac- 
tion Hp as 

N L/a -, ^ 

Hp ^ A^5|5-+, ^ A^(ata, - -)ia]^^a,^, - -) « 

L 

(47r^aA) ^ dxj {[: ^l^(a;j)t/'L(a;j) : + : ijli{xj)iljR{xj) : 



+e-^^'=--^V^i(a:,)V^«(a:,) + e'^''-^^i^l{x,)M^,)]x 
[: '0L(3^j+i)V'L(2^j+i) : + : '0i(^j+i)V'R(2^j+i) : + 

e-'''-^^^^i^[{x,+,)M^,+i) + e'''-^^^^i^U^,+,)M^j+i)]} = + H^p , (23) 
where 

L 

hP = (47r2aA) y" dxj [: '0i,(xj)V'L(xj) :: '0i,(xj-+i)V'L(xj-+i) : 



+ ■ ■^Lixj)ipL{xj) :: ilj\i{xj+i)iljR{xj+i) : + : ^]j(a;j)^R(a;j) :: ilj[{xj+i)ilJL{xj+i) :] , (24) 

and 

L 

H^? = (47r2aA) y do:,- [V^i(a:,)V^R(x,)e-2*^^^^ + 4(x,)V^L(x,)e2**=^^^] x 



8 



While evaluating Hp^ is rather straightforward, since it contains only normal- 

(2) 

ordered fermionic left- and right- densities, evaluating Hp is a little bit more 
involved, due to "crossed" L — i?-interaction. In fact, at any kp, momentum 
conservation selects the pertinent contribution to Eq.(25), given by 



-2ikFa„iA 



(26) 



where a possible "Umklapp" contribution, arising when /c^a ~ 7r/2, has been 
neglected ^ . To normal order Hp^ , one may rewrite it as 



(47r'aA)e'^'=^» / dx 



: 'ip[{x)ipL{x + a) : + 



27ra, 



27ra. 



(47r2aA)e-2^'=^"y dx 



27ra 



2'Ka 



,(27) 



which, for a ^ 0, becomes 



J-/ 

hP = -2(47r^aA) cos(2A;f a) j dx : ijj[{x)ipL{x) :: ij;\i{x)ijjR{x) : 



L 

+AnAsm{2kFa) j dx [: ipl{x)ipL{x) : + : '0lj(a^)V'fl(2^) :] 



dx ip\{x) — ^ i^^nix 



dx 



(28) 



^ This is the case, for instance, of the "half filled" fermionic sea in zero chemical 
potential 



The various terms in Eq.(28) may be interpreted as follows: 

• A shift in the chemical potential: 

L 

47rAsin(2/cFa) ^ dx [: ilj[{x)iljL{x) : + : :] , 



which is accounted for by simply redefining kp through the equation 
- (a J) cosikpa) + 2A sin(2A;i.a) = H ; (29) 

• A L — R interaction term: 

L 

— 2(47r^aA) cos(2A;i?a) j dx : tl)]^{x)tl)L{x) :: ■?/']j(a;)-?/'R(a;) : , 



that adds up to a similar term coming from Hp\ giving 

L 

2(47r2aA)[l - cos(2A;Fa)] J dx : ij[{x)ijL{x) :: ■ip\i{x)ijR{x) : ; (30) 



A renormalization of the Fermi velocity given by 

dx i^\{x) ; '0l?(^) 

dx 





dipR{x) 



dx 



(31) 



Using the well-known bosonization rules (A. 9), the fermionic Hamiltonian H^^ 
may be written in bosonic coordinates as 



H'^'-^^ldx 
An 



dx 



+ 



dx 



+ 2 



i± [ 



dx 



d(j)L d(l)R 
dx dx 



,(32) 



where g2 = 94 = 47r(aA)[l — cos{2kFa)]- 

One may readily see that H'' corresponds to the Hamiltonian for a free, mass- 
less, real bosonic field $ in 1-1-1 dimensions, which is described by the Hamil- 
tonian 



1^ 

H[U, [dx 

A.'K J 



g ^\dx 



(33) 



where the momentum conjugate to $ is H = (2t^/ g)^. 
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Upon defining two independent chiral fields, and 0^, as 



d(()'[{x ~ vt) 1 
dx ~ \f2 



dx 



1 
72 



27r ^ ^9$ 



(34) 



one immediately sees that 



H[IiM ^ H[4>IA'r] ^ ^ I dx 





1^ 

I' 



dx 



dx 



(35) 



which, when expressed in terms of 0£ ^ and 0^ ^, yields Eq.(32), provided 
that 



= \/ (vf + 92^ - gl ; 9 



I vf + 92 + Qa 
vf + 92 - 94: 



(36) 



Thus, the correlation functions of all the operators depending on (f)^ and 
may be evaluated by the replacements 

4>L-(t>R-^^gWL-<t>'R] , 0L + 0ii = + 0fe] , (37) 

with 0^ free, chiral bosonic fields. 



4 Phase diagram of the JJ-chain 



In Ref.[14], it has been evidenced that the phases allowed to a JJ-chain are: 

(1) A "Mott insulating" (MI) phase; 

(2) A "band isulating" (BI) phase; 

(3) A Repulsive Tomonaga-Luttinger phase (RTL); 

(4) A Superconducting (S), attractive Tomonaga-Luttinger phase. 

Here, we shall determine the range of the JJ-chain parameters associated to 
each allowed phase in the Vg — J-plane and, using the TLL approach, we shall 
provide a careful derivation of the phase boundaries; of course, our results 
crucially depend on the approximations made in Section 2. Our subsequent 
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analysis is based on the bosonic, low-energy effective Hamiltonian given in 
Eq.(32). 

To analyze the onset of the MI phase, one has to include also the Umklapp 
term Hp in Eq.(25), given by 

L 

^ -(aA) j dxj [i){{xj)i){{xj+i)^R{xj)i)R{xj+i) + ^^ji{xj)i)^R{xj+i)i)L{xj)ipL{xj+i)] , (38) 





whose bosonized version yields 

//^ = - (2|^) / dx : cos[2V2*(x)] : ; gu ^ «^(^)' " (3^) 





Eq.(39) and Eq.(33) yield the Hamiltonian for a 1+1- dimensional Sine-Gordon 
model, whose phase structure, as a function of the parameters g and gu has 
been extensively studied [19]. There are two distinct regimes: if < 2, the 
interaction is irrelevant and the theory is perturbative in gu-, while, ii g > 2, 
the interaction is relevant. In the thermodynamic limit (L oo) the system 
flows towards a strongly-coupled regime, where the Umklapp interaction is 
responsible for the creation of a gap in the excitation spectrum and for the 
onset of long range Ising-Neel order [7] . In the language of the J J-chain, this 
corresponds to a checkboard charge ordered state with the charge at each 
grain either equal to n, or to n -|- 1: this is the Ml-phase. 

The MI charge-ordered region in the Vg — J plane may be identified with the 
condition g > 2, which reads 

isin(ak,)U-lQ>lj . (40) 



As J = 0, Eq.(40) is satisfied by any value of kp (provided that there is a 
real solution to Eq.(29)) When there is no real solutions to Eq.(29), that is, 
for large enough \H\, the chain undergoes a phase transition towards the BI 
phase. This shows that, for J — 0, the only possible phases are either the BI 
phase, or the MI charge-ordered phase. 

To see how the transition towards the BI phase extends for J > 0, one may 
start again from Eq.(29), describing the booundary of the BI phase. If if > 0, 
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there are no real solutions of Eq.(29) for 

3 7^ 

H-2A>J =^ H -2Ez- J + -—>0 . (41) 

8 Ec 



As if < 0, on the other hand, there are no real solutions if 

3 P 

H + 2A<-J =^ H + 2Ez + J--—<0 . (42) 

8 Ec 



Eqs.(41,42) define two regions in the phase diagram corresponding to BI 
phases, since the density of charge-carrying states at the Fermi surface is 
0. 

Furthermore, as J > 0, Eq.(40) admits real solutions only if 



<1 ^ J <\ El + ^EzEc-Ec , (43) 



which implies that the MI phase closes when J = Ec + ^EzEc — Ec- Since 

A changes sign for J = J* = ^^EzEc, one finds that, as the MI phase closes, 
the Tomonaga-Luttinger liquid interaction is still repulsive (that is, (? > 1). 

The phase where, instead, the Tomonaga-Luttinger liquid is attractive (which 
is a necessary condition, to achieve superconducting correlations in the 1-d 
system) takes place for J > J*, that is, for gf < 1. 

In Fig.l we plot the phase diagram obtained using the TLL-approach. We 
observe that, due to the rcnormalization of the Fermi velocity, the line cor- 
responding to 5^ = 1 is a straight horizontal line: thus, as long as the TLL- 
description of the JJ-chain holds, one cannot push the system across this line 
by acting on the gate voltage Vg. We expect that this behavior is a byproduct 
of the approximations introduced in Section 2; higher order contributions to 
Hes should strongly modify the line corresponding to gi = 1. 



5 Two-boundeiry Sine- Gordon- mo del description of a finite JJ-chain 



In the following, we shall consider a one-dimensional JJ-chain with a weak 
link (i.e., a junction with a different nominal value of the Josephson coupling, 
Ew) at its center, whose position is set at a; = 0, and ending with two bulk 
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J 




Fig. 1. Sketch of the phase diagram of the JJ-chain in the Vg — J plane derived 
within TLL-approach, as discussed in Section IV. 

superconductors, whose phase difference is held fixed at (p (i.e., ipji = —<fL — 
(p/2). Using the bosonization method, in this Section we show that this finite 
JJ-chain is pertinently described by a two-boundary Sine-Gordon model [2]. 

Upon introducing JW fermions on both sides of the weak link, one gets 



where the labels > (<) refer to observables at the right (left)-hand side of the 
weak link. 

Using the long wavelength approximation, the fermionic string in the expo- 
nential of Eqs.(44), is easily evaluated as: 

L/a-l L-a 

1 = 1 Q_ 
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ttL 1 

— + -[0L,>(L-a) + 0R,>(L-a)] , (45) 



which, in turn, imphes: 



^L/a,> L J ■ ^ \2naJ 



1 



and, by keeping only the leading contributions to Eq.(46) in the cutoff a, as 
a — > 0, one gets, for x > 0, 



\27raJ ■ ■■ ■ ■ \ J 



Similarly, for x < 0, one obtains 

= e^['^^'<(-^)-'^«-<(-^)l . (48) 



The boundary condition at x = is, instead, dynamical, since it depends on 
the strength of the weak link, Ew In terms of the spin variables, the weak 
link interaction may be represented as a pointwise contact Hamiltonian given 

by 

^^^ = ^[^oVo:>+'5oVo:<] ■ (49) 



Taking into account that 5"^ = |[4ao + aoao] = f[- ^l(0)^l(0) +^/'i(0)V'i?(0) :] 
and the requirement that SqSq—SqSq — 2Sq, for any value of g, the operators 



and Sq are realized as: 



s 
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S-^^(^^y :e-^2[M0)-M0)], . (50) 



^From Eqs.(50), the dependence of Hw on the bosonic coordinates is given 
by: 



aEw /27ra 



w 



2L 



L 



exp 



;[-0i„<(O) + 0^,<(O) - 0z„>(O) + 0K,>(O) 



+h.( 
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Using Eq.(37), one immediately sees that the boundary interaction Hamilto- 
nian at x — takes the form 

= /27ray ^_ ^i^f [<^^,+(o)-0h,+(o)] . ^ . ^-i^ict>L,+{o)-'l>R,+m .](52) 

2^ \ L J 



where 



d(f)L,+{x - vt) 
dx 



.{x + vt) 



dx 



1 

1 

V2 



(53) 



with $±(x, = ^[$>(x, ± $<(-x, i)]. 

The boundary conditions at the bulk superconductors may be written as: 
^[ct>l,+ m t) - cl>%4L, t)] = <p (mod 2nk) . (54) 



By inspection of Eqs.(52), one sees that the field $_(a:, t) fully decouples from 
the weak link dynamics. Furthermore, its boTindary condition is $_(L, t) = 
Vt, thus implying that $_ is insensitive to variations in the phase difference 
between the bulk superconductors. As a result, the field $_ does not contribute 
to the dynamics of the JJ-network. Using Eq.(52), one gets that the pertinent 
effective Hamiltonian Hjj is given by 



H 



jj 



AttJ 



dx 



dx 



+ 



dx 



aE\v f2Txa\^ 

—) :cos 



:(55) 



with {ipL +, (f)R+) {(pL, <Pr), while the pertinent boundary condition is given 
by: 



^IML: t) - t)] = ip (mod 27rA;). 



(56) 



The model described by Hjj, supplemented with the boundary condition 
in Eq.(56), coincides with the two-boundary Sine-Gordon Hamiltonian intro- 
duced in Eq.(l), provided $ in Eq.(l) is identified with (0l — 0h)/-\/2, A^j is 
sent to oo, and A^, is identified with irr'Y ■ Intuitively speaking, while 
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the boundary condition at a; = L is always Dirichlet-like, ai x — 0, since Ew 
is finite, the boundary condition is dynamical, and is given by: 



aEw /2'KaY ^ 



|(</'4o,t)-Mo,t)) 



+ ^;^[0L(O,t)-0ij(O,t)]=O .(57) 



In particular, for small Et^^ Eq.(57) provides Neumann-like boundary condi- 
tions for 01, — at a; = 0, while it provides Dirichlet-like boundary conditions 
for large values of Ew- 

For gi < 1, the boundary term is a relevant operator and one may use the 
Renormalizion Group to describe how the renormalization of affects the 
ground-state energy: as we shall evidence in section 6, there is a renormal- 
ization group invariant length scale L* such that, for L > L* the JJ-chain 
behaves nonperturbatively in Ey/ while, for L < L*, it behaves perturba- 
tively. At variance, for > 1, the JJ-chain is always perturbative in E^, since 
the boundary term is now an irrelevant operator. For g = 1, the boundary 
term is marginal and the bulk system is fully described by a pair of nonin- 
teracting chiral fermions and the partition function may be computed exactly 
[2]. 

To summarize, for > 1 E^ always flows to 0, while, for < 1, there is a 
characteristic " healing" length L* — L (^ Ew{^=i) ) ' separating a small- E'^t^ 
perturbative regime from the nonperturbative one Ew > 1- Similar features 
are exhibited by a superconducting loop closed by a Josephson junction of 
strength Ej, when Ej is regarded as an effective coupling strength [20]. 

In the next Section, we shall prove that the behavior of the DC Josephson 
current as a function of (p depends crucially on whether E\y flows to zero, or 
to large values. Namely, we shall flnd that, when Eyy flows to zero, the DC 
Josephson current has a sinusoidal behavior, while, when Ew ~ 1, one gets 
the sawtooth behavior. 



6 Josephon current across the JJ-network with a weaik link 



In this Section, we shall compute the functional dependence of the DC Joseph- 
son current on the phase difference between the bulk superconductors, by eval- 
uating the zero temperature canonical partition function Z[Ew], from which 
the Josephson current may be evaluated as 

c o(p 
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with 



Ejjy\ = - lim i In 



Z[E- 



w\ 



(59) 



The partition function of the two-boundary Sine-Gordon model has been ex- 
actly computed for particular values of g [2] ; due to our interest in providing an 
estimate of the Josephson current fon any value of g^, we resort to an approxi- 
mate computation based on the Coulomb Gas Renormalization Group scheme 
[21], described in detail in Appendix B. Our analysis shows that, while for Neu- 
mann boundary conditions {Ew 0), one gets I{ip) oc sin((^), for Dirichlet 
boundary conditions {Ew > 1) one gets a sawtooth dependence of on ip. 



Case A: Ew 



As Ew ~ 0, Eq.(57) provides Neumann-like boundary condition at x = 0, 



(60) 



together with the Dirichlet-like boundary condition x — L: 



(f)L{L-vt)-(l)R{L + vt) 



Using the mode expansion given in Eq.(A.5), one sees that Eq.(60) is satisfied 
if 

f _ 
Ql-Qr^— ; pl^Pr = P 

V9 



(.L(n)e^'="^ = -<f>R{-n)e-"'-'' = a{n) ; (fc, = ^ + ^) , n e Z) (61) 



which, in turn, imphes 



(o 47r 

01,(0 - vt) - + = ^ + 1- E 



-iknVt 



-a{n) = x{t) . 



(62) 
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The relevant vertex operators are given by 



with 



p-iknVt p-iknVt 

x+W = ^ E ^— «(^) ; x-{t) = - E ^— «(^) ■ (64) 

^ n<0 '^"^ n>0 



The partition function is given by 



cxp 



aEw f2na\^ 



y o?T : cos[^|x(t)] : 



(65) 



where t — it, Tt- is the (imaginary) time-ordered product, and the brackets 
(. . .) mean that expectation values should be computed with respect to the 
ground state of the Hamiltonian in Eq.(35). 



Using the expansion of Eq.(65) in a power series of Ew as 

Z[E^] = Z[E^ = 0] f: E Jtldn (tAU : e^^'^'^^^'^^ :]) (66) 



a,=±l n k=l \ k=l 



with Ew^^ {^tY^ and using the identity 



^[ri-r2][x-(ri),x+(r2)]=41n[l-e-^"^^] , (67) 



one gets 

, - s ■ P Ti-a/v Tj-i-a/v 







Z[Ew = 0] 



i=0 ^ a3=±l 

2i 
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since Wick's theorem implies ^ 



2j 

fe=i 



2i 

n [i-e^ 

u<r=l 



(69) 



To analyze the short-distance divergences in Eq.(68), one has to rescale the 
cutoff a a/K, A > 1, in order to approximate nu<r=i[l ~ e~E^I^"~^'"l]^^"""'' 



with n 



2j 

u<r=l 



2gauOCr 



. As a result, one obtains 



/3 Ti—a/Av ij — \—u,/i\u 2j 

fdn I dT2... I (ir,e'^^'^=i"'= n 

n n n u<r=l 



2tiv 



2ga„ar 



(3 Ti-a/v (J 2j 

on n u<r=l 



Tj — i—a/ V 



2ttv 



2gauCtr 



At a given order 2j, and for A — > oo, the most diverging contributions come 
from integrals containing an equal number of positive and negative a's. Thus, 
A-scaling of the integrals appearing in Eq.(70) is taken into account by means 
of a multiplicative renormalization of the effective coupling strength 

Ew Ew{A-) = A^-^Ew{^ = 1) . (71) 



Eq.(71) implies that the boundary interaction at x = is irrelevant and 
the Neumann fixed point is always stable for g > 1 (i.e., in the repulsive 
Tomonaga-Luttinger phase). The RTL phase is always associated to a stable 
Neumann fixed point. 

To evaluate /((/?), one may retain only the first order terms in the Ew 
expansion in Eq.(66), getting 

- P 



^ In Eq.(68), the cutoff a has been introduced to regularize possible short-distance 
divergencies in the argument of the integral. It should be identified with the lattice 
step introduced in Eq.(2). 
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from which the network energy is derived as 

£,,M = -ltail.^ = (a^,)(2p)"co.(^) . (73) 

Using Eq.(58) one gets /(</?) oc sin((/7). 

It is comforting to see that Eq.(73) reproduces the pertinent renormahzation 
of the effective couphng constant given in [14] . 



Case B: Ew 1 



When analyzing the case in which the effective couphng grows, as the size 
of the system goes large, wc must consider that our analysis started from 
expanding fermionic fields whose band energy is equal to J. Accordingly, the 
scahng should stop as Ew ~ J- The scale at which this happens, A*, is found 
by the condition 

Ew{A=l){A*y'' = J (74) 



This implies that scaling stops as the size of the system becomes of order of 
L*, given by 



For L < L*, the theory is still perturbative. As L ~ L*, instead, the system 
enters the nonperturbative region. In this limit, the field 0l — 0ij has to satisfy 
Dirichlet-like boundary conditions both at a; = L and at x — 0. Prom Eq.(57), 
one gets 

sin[^|[0i(O - vt) - 0^(0 + vt)]] = -_— - — [0^(0 - vt) - 0^(0 + vt)] 

thus, the Dirichlet-like boundary condition at a; = is 

0l(O - ^;^) - 0ii(O + t;^) = . (77) 

Using the mode expansion given in Eq.(A.5), Eq.(77) may be easily satisfied 
by setting 

qL = qR = q ; Pl=Pr = P 
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ik„L 



(78) 



provided that: 



27rn 



kn — ; n E Z ; —^^-Kp = 2nk + (p 
Ij 



(t>L{n) = -(pni-n) = a{n) . 



(79) 



The partition function is now given by 
Z = Trexp 



n>0 



with the pertinent boundary conditions given by Eq.(79). 

The trace in Eq.(80) may be factorized into a contribution from the oscillatory 
modes, and a contribution from the zero modes, so that 



Z Zaar Zo- 



'osc-^0— modes 



(81) 



with = [n~ i[l - [e-^'^r]]-', and 



oo 



-2o-modes — 2^ CXp 

A;=— oo 



'^2gL 2ttJ 



(82) 



^From Eqs.(81,82), one gets /((/?) oc (p, for — tt < (p < tt; this yields the 
well-known sawtooth behavior [20]. 

The switch to this behavior from the sinusoidal behavior obtained for Ew ~ 
signals the crossover from a perturbative to a nonperturbative regime of the 
chain. It should be observed that, for g > 1, the DC Josephson current has 
always a sinusoidal dependence on the phase difference between the bulk su- 
perconductors since, in this region, the boundary term is an irrelevant operator 
and, thus, the chain's behavior is always bulk-dominated. 
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7 Concluding remEirks 



Our analysis shows how a two-boundary Sine-Gordon model emerges as a per- 
tinent effective description of a finite JJ-chain. Since our analysys heavily relies 
on the bosonization method, we expect that boundary Sine-Gordon models 
may be useful where the Tomonaga-Luttinger liquid paradigm is relevant. For 
instance, a magnetic spin system with a pertinent impurity at its center and 
with the spins at its extrema held fixed, may support a spin current across the 
chain, with different behaviors, depending on the boundary conditions around 
the impurity. Similarly, one may envisage other applications of our results to 
quantum wires [22] and carbon nanotubes [23]. 

According to the g^-theorem [13], the boundary entropy of the chain should 
always decrease, as one gets towards the thermodynamic limit. Thus, for L > 
L* — > oo, one has 

Sd - Sn ^ lim [Urn [In Zd -In Zn]] , (83) 

/3— >oo L— >oo 



where Zd/n is the partition function computed with Dirichlet / Neumann 
boundary conditions, respectively. Eq.(83) yields a nonvanishing result only 
because of the contribution of the zero modes; namely, from Eq.(82), one gets 
that 

SD-Sr,^H^^ ■ (84) 



Remarkably, the entropy variation depends on the sign of In y^. Thus, as 
g > 1 (i.e., within the RTLL phase) the Dirichlet boundary entropy is higher 
than the Neumann boundary entropy and then, in the thermodynamic limit 
the system flows from the Dirichlet to the Neumann fixed point. Gonversely, 
as g < 1 (i.e., within the superconducting region) the Neumann fixed point 
carries an entropy that is higher than the one associated to the Dirichlet fixed 
point. Accordingly, the flow now goes from the Neumann to the Dirichlet fixed 
point. 

The evaluation of the Josephson current from the partition function of the 
two-boundary Sine-Gordon model explicitly shows, for g < 1, a crossover 
from a perturbative regime {E^y ~ 0), in which the current is a sinusoidal 
function of the phase difference at the boundary, ip, to a nonperturbative 
regime {Ew/J > 1), where it exhibits a sawtooth functional dependence on 
if. 

There is a striking, and yet intuitive, similarity between the finite JJ-chain 
investigated in this paper, and an rf-SQUID in an external magnetic field [20]. 
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For the latter system, a variation of the flux threaded by a superconduct- 
ing loop operates the crossover in the behavior of the Josephson current, as 
a function of the applied flux. This might suggest that other very interest- 
ing condensed matter realization of boundary Sine-Gordon models, may be 
provided by superconducting loops, interrupted by two, or more, Josephson 
junctions [24]. Remarkably, using the results in Section 6, one immediately sees 
that the effective potential of the finite JJ-chain, as a function of the phase 
difference at the boundary of the weak link, i.e., ip = (0l(O) — 0r(O)), exhibits 
only one minimum within the RTLL phase (i.e. ijj = (p), since ~ 0, while 
it is a two- level quantum system in the superconducting region [Ew/ J > 1), 
provided that (/? ~ tt. 
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A Bosonization rules 



In this Appendix the bosonization rules used in this paper are reviewed. 

It is a peculiar property of 1 + 1 dimensional theories that it is possible to 
realize chiral fermionic fields in terms of chiral bosonic fields, and vice versa. 
If one starts, for instance, from the chiral components of a free, massless, 
Klein-Gordon field $ in 1-1-1 dimensions, the equation of motion for $ is 

$(x,i) = , (A.l) 



92 2 
V 



where periodic boundary conditions are assumed. $ may be written, as the 
sum of two chiral flelds as 

$(x, t) = -^['Pl{x - vt) + + vt)] (A.2) 

with (f)L and chiral Fubini-Veneziano flelds [25], whose mode-expansion is 
given by 

27r 2m ^iknix-vt) 
(pL{x,t) = (f)L{x - vt) ^Ql- -rPL{x - vpt) + —j-Yl Z 0l(?T')(A.3) 
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and: 

fhr,(nr. f.\ = rhr,( nr. A- iif.\ = rir, -I rtrtf t. A- iit\ 4- - 

L ^ k 



2n 2ni ^iKix+vt) 
r{x, t) = (f)R{x + vt) = qR+ —Pr{x + vt) + — ^ (j)R{n) .(A.4) 



The basic commutation rules are 

[Ql:Pl] = [qr^Pr] = i ; [0L(n),0L(m)] = n5n+m,o ; [0RN,0i?M] = -n5n+m,o , (A.5) 

and {kn} is a (discrete) set of nonzero modes depending on the boundary 
conditions imposed on the bosonic fields (for instance, for periodic boundary 
conditions, one gets fc„ = n e Z). 

Due to the commutation rules in Eq.(A.5), the bosonic vacuum |Bos) is defined 
by 

Pi I Bos) ^pr\Bos) = 

Mn)\Bos) = 0R(-n)|Bos) = (n > 0) , (A.6) 

and one may then define a creation and an annihilation part for each field, 
i.e., 

(t>L{x) ^qL + -j- Yl -[—^L{n) ; (j)l{x) = — -xpL + -j- -r—<l>L{n){A.7) 



and 



^ kn>0 '^^ 



<t>R{n) ; (I)r{x) ^^xpR+'^ Y '-j—<t>R{n) .(A.S) 

^ ^ k„<Q 



^Prom the commutators given in Eqs.(A.5), one sees that the modes of the 
operator obey the same algebra as the modes of the fcrmionic density 

operator, normal ordered with respect to |FS). Thus, the fermionic bilinear 
density operator : ipiipL '■ may be identified with the bosonic density operator, 
27^^"' P'^ovided that |Bos) is identified with |FS). The same identification may 
be carried for the ii!-modes. Therefore, one gets a first bosonization rule 
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The second rule is obtained if one identifies the chiral fermionic fields with 
normal ordered vertex operators of bosonic fields, defined by 



V27ra 



(A.IO) 



The correspondence rules are now given by 



i)i{x - vt) 



: e" 



; tPU^ + vt)^^: e-^'^«(^+''*) : .(A.ll) 



To check the consistency of Eq.(A.ll), one has to consider the "braiding rule" : 



: e 



i<x<t>L/R{x) .. „il3(t>L/R{y) .— J-^a/B . „iP(l>L/R{y) .. pia'l>L/R{x) . 



(A.12) 



and the vertex- vertex correlators: 



(Bos| : e^'*^^/«(^T''*) :: g-^^'^^/s^^'T''*') : |Bos) 



±i 



2 J sm[^{x - x' ^ v{t - t'))] 



^Prom Eqs.(A.12,A.13), one derives the basic anticommutators: 

— vt),ip\^{x —vt)} = 5[x — x —v{t — t)] 

and 



{^R{x + vt),tlj]i{x' +vt)} = 5[x - X + v{t-t)] 



(A.14) 



The other correlators used in the paper are derived from Eq.(A.13) and from 
Wick's theorem applied to normal ordered vertices [25]. 



B Renormalization Group equations for the JJ-chain with a weeik 
link 



In this Appendix, the flow of E\y is derived within the Coulomb Gas renor- 
malizion Group approach. 
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To analyze the Renormalization Group flow for the JJ-chain with a weak 
link, one needs to observe that, at order 2j, the short-distance most diverging 
contribution to the partition function is given by: 



Ti-a/v T2j-i-a/v 



dT2j 



2j 

S n [1 - e"^!^""^"!]^^"""'' (B.l) 

ai+...+a2j=0 u<r=l 



Rescaling the short-distance cutoff as a — > a/A and sending A — > oo implies 
the following renormalization group scaling equation for Ew: 

M=(l-.)^.(A). (B.2) 



1 

1-9 



/,From Eq.(B.2), one may easily identify the cutoff scale A* = [ ^^^(^=1) 
at which Ew becomes ~ J. It means that a system of size L, with a weak link 
of nominal strength E^y = E]y{A = 1), "crosses over" towards the strongly 
coupled regime, as its size is increased to L* — h*L [14]. 

In addition, it has to be noticed that, to remove the cutoff, one needs a further 
renormalization, due to "one-dimensional charge annihilation processes" . This 
may be evidenced, for instance, by applying Anderson- Yuval-Hamann analysis 
of a one-dimensional instanton gas [21]. 

As the cutoff is rescaled from a to a/A, two charges, of opposite sign, may 
annihilate with each other, if they were originally separated by a distance 
between a/{vA) and a/v. As a result, the integral at order 2j + 2 should 
provide an extra contribution to the integral at order 2j, which we are now 
going to calculate. 

Upon defining 

r = ^±^ ; r = r+-r_ (B.3) 



where t+ being the coordinate of the -|-l-charge, and r_ the coordinate of the 
— 1-charge, the extra contribution arising to order 2j, is given by 



a/v T\—a/L T'zj — i^a/v 72j izj—i—u^rv 

(Ewf^^^ J dT^jdn j dT2... J dT2jX JdT+ J dT + 

a/iAv) 



T2j r2j-i-a/v 

■ / 

T2j-2-a/v 
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2nVT 



2g 



2j 

n [l-e-^l""-"'-l]25"""'-exp 

u<r=l 



a/v f) Ti—a/L T2j-i—a/v 

+{Ewf'^^ I dT^jdn j dT2... I dT2, 
a/{Av) 



T2j-i-a/v T2j-2-a./v 

dT+ J dT + . 

T2j-a/v T2j-3-a/v 



I ^^^^ I I 



1 



2i 



2-KVT 
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n [1 

w<r=l 



— p-^|Tu-TT-|l2sra„Or 



exp 



(B.4) 



If one expands the exponentials as 



exp 



k=l 



k=l 



one may derive the renormahzation group equation for ^ . The result is [21]: 

L 



g^g + dg^g- 



2TTva 



g[Ew] dlnA 



(B.6) 



As expected [19], the wavefunction renormahzation is needed only for g < 1. 

This completes the renormahzation scheme derived within the perturbative 
approach for the system in the presence of a weak link. 
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